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Abstract 

We present numerical exact results for the polaronic band structure of the 
Holstein molecular crystal model in one and two dimensions. The use of direct 
Lanczos diagonalization technique, preserving the full dynamics and quantum 
nature of phonons, allows us to analyze in detail the renormalization of both 
quasiparticle bandwidth and dispersion by the electron-phonon interaction. 
For the two-dimensional case some of our exact data are compared with the 
results obtained in the framework of a recently developed finite cluster strong- 
coupling perturbation theory. 
PACS number(s): 71.38.+i 
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The very fundamental problem of a single conduction electron coupled to bosonic degrees of 
freedom is still not completely understood. In the case where the bosons are lattice phonons 
that is what is usually called the polaron problem and numerous analytical techniques were 
tried out to tackle this complicated many-body problem in terms of the simplest electron- 
phonon (EP) Hamiltonians ]]J. Depending on the relative importance of the long-range or 
short-range electron-lattice interactions simplified models of the Frohlich g] or Holstein || 
type, respectively, have been widely used to analyze polaronic effects in solids with displace- 
able atoms. The crucial point is that the usual phase transition concept fails to describe the 
crossover from an only weakly dressed electron to the strongly renormalized less-mobile po- 
laronic quasiparticle with increasing EP coupling strength. Furthermore, as yet none of the 
various analytical treatments, based on variational approaches 0] or on weak-coupling || 
and strong-coupling adiabatic || and non-adiabatic |||7j perturbation expansions, are suit- 
able for the investigation of the physically most interesting transition region, where the 
highly non-linear "self-trapping" process of the charge carrier takes place. That is because 
precisely in this situation the characteristic electronic and phononic energy scales are not well 
separated and non-adiabatic effects become increasingly important, implying a breakdown 
of the standard Migdal approximation ||. In principle, quasi approximation-free numerical 
methods like (quantum) Monte Carlo simulations || and exact diagonalizations (ED) fH| 
can close the gap between the weak and strong EP coupling limits and therefore, at the mo- 
ment, provide the only reliable tool for studying polarons close to the crossover regime. Along 
this line, previous ED work has concentrated on spectral properties of the one-dimensional 
(ID) Holstein polaron. In practice, however, memory limitations have restricted the calcu- 
lations to either small values of the EP coupling constant or to very small clusters fll|-|13|, 
thus, e.g., making a more detailed discussion of the polaronic band dispersion impossible. 
First ED studies performed on the ID ten-site Holstein model indicate that the small po- 
laron dispersion may differ significantly from a rescaled cosine tight-binding band especially 
in the crossover regime []14|| ; a result which has been corroborated more recently by a finite 
cluster (FC) strong-coupling perturbation theory (SCPT) fl5| . 



Encouraged by these findings, it is the aim of the present paper to perform a systematic 
Lanczos study of the ID and 2D Holstein model on large enough lattices, in order to discuss 
the evolution of the polaronic quasiparticle band structure in dependence on the phonon 
frequency and EP interaction strength. 

The Holstein Hamiltonian reads: 

h = -tJ2( c h + c \ c ^ - v / %^E( fe l + b i)cUi + + 1) > (!) 

(ij) * » 

where cf and bf are the annihilation [creation] operators of a (spinless) fermion and a boson 
(phonon) at Wannier site i, respectively. In (1), the following idealizations of real electron- 
phonon systems are made: (i) the electron transfer t is restricted to nearest-neighbour (NN) 
pairs (ij); (ii) the charge carrier is locally coupled to a dispersionsless optical phonon mode 
(e p denotes the EP coupling constant and hu is the bare phonon frequency); (iii) the phonons 
are treated within an harmonic approximation. Recall that the physics of the Holstein model 
is governed by two dimensionless parameters A = e p /2Dt and a = Je p /?kv; i.e., irrespective 
of the adiabaticity ratio 7 = t/hu, small polarons are formed provided that A > 1 and a > 1. 

In the numerical analysis of the Holstein model, we employ the standard Lanczos al- 
gorithm in combination with a well-controlled truncation of the phononic Hilbert space. 
According to the truncation procedure described in [f[4| , we have to check for the conver- 
gence of both the ground-state energy E (M) and the phonon distribution function, which 
has to be determined self-consistently in the ground state \\Pq(M)). Increasing the maximum 
number of phonons retained (M), convergence is assumed to be achieved if the relative error 
\Eq{M + 1) - E {M)\/\E {M)\ becomes less than 1(T 7 . 

In order to discuss the single-polaron polaron band structure, we make use of transla- 
tional invariance and compute the wave-vector resolved spectral density function 

AAE) = J2\(^\c^\0)\ 2 6(E-E^), (2) 

n 

using a polynomial moment expansion together with the maximum entropy method | |16| . 
Then the so-called "coherent" band dispersion, Eg, is obtained from the first peak of Ag(E) 
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having finite spectral weight 

Z k = K^oJl c i?l°)| 2 ' ( 3 ) 

where ^) denotes the single-electron state being lowest in energy in a certain if-sector. 

In the first place, we focus on the ID model and start with the study of the weak-coupling 
case (A < 1). Figure 1 displays the band dispersion calculated at e p = 0.1 for various 
phonon frequencies ranging from the adiabatic (7 = 10) to the non-adiabatic (7 = 0.125) 
regime (in the following all energies are measured in units of t). For low and intermediate 
phonon frequencies, the energy to excite one phonon lies inside the bare tight-binding band, 
l cos Ki } of a D-dimensional hypercubic lattice with unit lattice spacing. Thus 
at arbitrarily small e p , predominantly "phononic" states enter the low-energy spectrum in 
those ^-sectors for which Eg being separated from the ground state by an energy E ~ 
hu. As a result the dispersion curves reflect the well-known peculiarities of the absorption 
spectra of an electronic or excitonic system weakly interacting with dispersionsless optical 
phonons P!7|JI^| , i.e., a nearly unaffected cosine dispersion near the band center (K = 0) 
and a practically flat region at larger momenta. Note that even in cases where the phonon 
frequency is comparable to the bare electronic bandwidth a remarkable "flattening" of the 
band structure appears in the vicinity of the zone boundary (e.g., at ftiu — 4 we found 
E n (e p = 0.1)/i? 7r (0) ~ 0.92). In the extreme non-adiabatic regime (7 <C 1) the phonon 
distribution function is sharply peaked at the zero-phonon state and, because the interaction 
is weak and the phonons can follow the electron instantaneously, we obtain a negligible 
renormalization of the band structure. To substantiate this interpretation we have calculated 
the .fT-dependent mean phonon number, 

i 

in the lowest band states l^ 1 ^)- The results are presented in the inset of Fig. 1 for three 
characteristic momenta K = 0, OAir, and n. Evidently the states keeping small momenta 
are basically zero-phonon states. For larger momentum, we observe a gradual crossover 
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from one-phonon to zero-phonon states upon increasing the phonon frequency, whereby the 
admixture of "electronic character" increases. 

Of course it is of special interest to understand the evolution of the band structure by 
increasing the EP coupling strength. Thus it is fruitful to compare Eg for the weak and 
strong-coupling cases at fixed phonon frequency. This has been done in Fig. 2 for various 
system sizes N = 10, . . . ,20, using EP couplings e p = 0.5 (a) and e p = 4.0 (b) together 
with an intermediate adiabaticity parameter 7 = 1.25. First of all one notices that the 
finite-size effects are very small, i.e., although the data points belong to different lattice 
sizes we obtain a remarkably smooth behaviour of Eg. As might be expected, for e p = 0.5 
the coherent bandwidth, AE = supgEg — infgEg, is approximately given by the phonon 
energy (AE = 0.782 ~ hu = 0.8) and the dispersion again shows the coexistence of two 
different types of band states at small and large momenta. By contrast, if the EP interaction 
is enhanced, small polaron formation takes place and we observe all signs of the famous 
polaronic band collapse. However, note that even in the relatively strong EP coupling regime 
displayed Fig. 2 (b) the standard Lang-Firsov formula, AE^p = 4Dt exp[— e p /huj] (obtained 
by performing the Lang-Firsov canonical transformation || and taking the expectation value 
of the kinetic energy over the transformed phonon vacuum), does not give a satisfactory 
estimate of the bandwidth. So we found AElf = 0.0269 which has to be contrasted with 
the exact result AE = 0.0579. Besides the band narrowing effect, there are several other 
features worth mentioning for polaronic band states in the crossover region. Most notably, 
throughout the whole Brillouin zone the band structure differs significantly from that of a 
rescaled tight-binding (cosine) band containing only NN hopping. As can be seen from a 
least-squares fit to an effective band dispersion Eg = J2f=o a i coslK, the residual polaron- 



phonon interaction generates longer-range hopping terms fl^ , [19|| . Concomitantly, the mass 
enhancement due to the EP interaction is weakened at the band minimum. It is important 
to realize that these effects are most pronounced at intermediate EP couplings and phonon 
frequencies. In this parameter region even higher-order SCPT, with its internal states 
containing some excited phonons, seems to be not tractable because the convergence of the 



series expansion is poor ||15|| . The deviation from the cosine dispersion at different e p is 
depicted in the inset of Fig. 2 (a), where we have used the eight and ten site lattices in 
order to reach larger EP coupling strengths. Here Eg is scaled with respect to the coherent 
bandwidth which strongly depends on e p and hu, for example, for N = 8 and hu = 0.8 
we found AE(e p ) = 0.1957, 0.0559, and 0.0143 for e p = 3.0, 4.0, and 5.0, respectively. 
It is obvious and doesn't need discussion that the dispersion is barely changed from a 
rescaled tight-binding band in the very extreme small polaron limit (A > 1). To visualize 
in more detail the different nature of large and small polaron band states, we have evaluated 
numerically the quasiparticle residue Zg and the mean phonon number Ng h ^ at all K- 
points. Obviously, the ./^-dependent renormalization factor can be taken as a measure of the 
"electronic contribution" to the polaronic quasiparticle. The results are shown in Fig. 2 (c). 
For weak EP coupling, one recognizes immediately the electronic and phononic character of 
the states at small and large momentum, respectively. The crossover between these two types 
of states is accompanied by a change in the mean phonon number by approximately one. 
With increasing e p a strong mixing of electron and phonon degrees of freedom takes place, 
whereby, forming a small polaron, both quantum objects completely loose their own identity. 
As expected, this leads to a significant suppression of the quasiparticle residue for all K. At 
large e p , the polaronic state is characterized by strong on-site electron-phonon correlations 
making the quasiparticle susceptible to self-trapping. Thus the coherent movement of the 
carrier is greatly hindered but itinerant band states do still exist. The small polaron is heavy 
because it has to drag with it a large number of phonons in the so-called "phonon cloud" 
(cf. Fig. 2 (c) inset). Since the high-energy band edge states are more vulnerable to decay 
from EP interaction the particle loses its spectral weight more rapidly at K = tt [^0] . 



Next, we consider the evolution of the band structure for the two-dimensional case. 
Figure 3 presents the results obtained for Eg and Zg along the highly symmetric directions 
of the Brillouin zone in the weak (a) and intermediate (b) coupling regime. As in the (weak- 
coupling) ID case, we found a flattening of the band structure in regions where E^± becomes 
comparable to the phonon energy. The sharp drop in the quasiparticle weight obtained with 



increasing \K\ at small huo is weakened in the non-adiabatic regime. That is, the transition 
region (0.25 ~ Z# ~ 0.75) is enlarged. In order to make a comparison with recent predictions 
of both second-order Rayleigh-Schrodinger SCPT and the re-summed FC-SCPT developed 



by Stephan [fig] , we have plotted Eg and Zg for e p = 5.0 and huj = 1.5 in the lower part of 
Fig. 3 . Interestingly enough, the flattening of the dispersion at the edges of the Brillouin 
zone survives to surprisingly strong EP couplings (note that Tiu ~ 2AE). The wave- 
vector renormalization factor Zg clearly indicates the polaronic nature of the charge carrier. 
That means, calculating the strong-coupling wave-function we have to include multiphonon 
states. Once again the Lang-Firsov formula (first-order SCPT) underestimates the exact 
bandwidth by more than a factor of two; AE^p = 0.2854. Even if we include second-order 
corrections the SCPT band structure significantly deviates from the exact result. Obviously, 
the agreement is poorest at the T-point (K = 0). Thus we are lead to the conclusion that 
the standard strong-coupling approximation fails to accurately predict the ground state 
energy and the polaronic bandwidth in the moderate strong-coupling regime. By contrast, 
the FC-SCPT approach fi~5|l , which in our case makes use of information extracted from a 
sequence of finite-cluster diagonalizations up to five sites, provides a correct description of 
the coherent band structure. The extremely good success of the FC-SCPT at intermediate 
EP couplings and phonon frequencies may be attributed to the inclusion of longer-range 
correlations, which are especially important in the crossover region from nearly free to small 
polarons. 

To summarize, we have performed a comprehensive study of the ID and 2D Holstein 
model using exact diagonalizations. Concerning the polaronic band structure of the Holstein 
model, our main results are the following. 

(i) In the weak-coupling case and for phonon frequencies less than the bare electronic 
bandwidth, we recover the expected flattening of the band dispersion at large momentum. 
Here our numerical results for the wave-vector resolved quasiparticle weight factor and the 
mean phonon number clearly demonstrate the electronic and phononic character of the band 
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states near the band center and band edges, respectively. 

(ii) Increasing the electron-phonon coupling the flattening survives to a large extent the 
crossover from large to small polarons. The small polaron state is characterized by less 
spectral weight and a large average number of phonons. 

(hi) At intermediate EP coupling strengths and phonon frequencies the effective pola- 
ronic band dispersion deviates substantially from a simple tight-binding cosine band due to 
further than nearest-neighbour ranged hopping processes generated by the residual polaron- 
phonon interaction, implying the importance of multiphonon states. In the intermediate 
coupling regime, the polaronic mass enhancement can be one order of magnitude smaller 
than predicted by standard strong-coupling perturbation theory. 

(iv) In the extreme strong-coupling limit (e p ^> ttiv, t), we have a nearly dispersionsless 
small polaron band with almost vanishing spectral weight. That means, here the exponential 
band narrowing completely dominates the physics of the Holstein model and, in real solids, 
coherent polaron motion will be immediately destroyed by disorder effects. 

The authors are greatly indebted to W. Stephan for putting his 2D FC-SCPT data at 
our disposal. H. F. is especially grateful to D. Ihle and J. Loos for many useful discussions. 
This work was performed under the auspices of Deutsche Forschungsgemeinschaft, SFB 
279, Bayreuth. Special thanks go to the LRZ Munchen for the generous granting of their 
computer facilities. 
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FIGURES 



FIG. 1. Polaronic band dispersion Eg for the ID Holstein model. Data points are exact 
results obtained for a finite lattice with N = 20 sites and at most M — 10 phonons using 
periodic boundary conditions. Here and in the following the energy scales of the ordinates 
were shifted by Nhu/2. The inset shows the mean phonon number Ng in band states with 
wave-vector K. 

FIG. 2. Polaronic band structure of the ID Holstein model in the weak (a) and strong (b) 
coupling regimes, where the phonon frequency is Two = 0.8. For both cases, the quasiparticle 
residue and mean phonon number are given in (c). 

FIG. 3. Band dispersion Eg (filled symbols) and quasiparticle weight factor Zg (open 
symbols) for the 2D Holstein model with e p = 0.5, Hu = 0.8 (circles), 4.0 (squares) and 
e p = 5.0, Tuo = 1.5 (diamonds). ED results are obtained using finite square lattices with 
N = 16 and 18 sites (the lines are only a guide for the eye). In the weak-coupling case (a) 
the chain dashed curve gives the free electron dispersion E^; for the strong-coupling case 
(b) exact data (filled diamonds) are compared with the predictions of standard second-order 
SCPT (dashed line) and FC-SCPT (solid line) pf. 
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Figure 1: 
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